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1. Introduction 

Various transport coefficients associated with a particle moving in low-dimensional 
periodic array of scatterers are often very irregular functions of control parameters. 
For example, the diffusion coefficient of a particle moving in the flower-shaped billiard 
(a 2D Lorentz gas with periodically distributed scatterers of six-petal flower shape) was 
found to be a very irregular function of the petal curvature [1] , and the conductivity 
of the classical Lorentz gas (with disc scatterers) varies rapidly with the applied field 
[2] . The reaction rate and the diffusion coefficient in the reaction-diffusion multibaker 
model also turned out to be irregular, fractal-like functions of the map's shift [3]. 
Similar behaviour was found for the diffusion coefficient of a particle thrown onto a 
periodically corrugated floor and subject to various types of external force, including 
Hamiltonian [4] and dissipative [5, 6] systems. In chaotic systems with two degrees 
of freedom the modes governing the relaxation to the thermodynamic equilibrium 
form a fractal structure with a nontrivial fractal dimension which can be related to 
the Lapunov exponent and the diffusion coefficient of the system [7, 8, 9]. Another 
example is provided by a deterministic dynamical system with dynamics defined by 
iteration of a one-dimensional map: the mean velocity and the drift coefficient in 
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this system were shown to be irregular, nowhere differentiable functions of the system 
control parameters [10, 11, 12, 13, 14]. 

It is natural to ask what the origin of these "irregularities" is, how to describe 
them quantitatively, and to what extent they are universal. Answering all these 
questions is far from trivial. One of the reasons why irregular behaviour of transport 
coefficients is so difficult to investigate is that until very recently our knowledge 
about this phenomenon was based mainly on numerical simulations. Although several 
different techniques of calculating these quantities were developed, e.g. the transition 
matrix technique combined with the escape-rate formalism [6, 13, 14, 15, 16, 17, 18], 
the Green-Kubo formula [6, 14, 15], or the periodic-orbit formalism [15, 19], they all 
lead to complicated and time-consuming numerical calculations of limited accuracy. 
Moreover, quite often these techniques are applicable only for some special values of 
the control parameters. What is worse, as a model becomes more realistic, the results 
are getting less accurate. For this reason the most precise data have been obtained 
for very simple, one-dimensional "toy models", while the question whether a more 
realistic 2D Lorentz gas (with scatterers of disc shape) exhibits any irregularities in 
the transport coefficient as functions of the disc radius is still a subject of controversy 
[6, 20]. 

This situation improved considerably when Groeneveld and Klages [10] gave 
exact formulae for the transport coefficients in a simple onc-dimcnsional deterministic 
dynamical system, introduced by Grossmann and Fujisaka [21], and defined by a piece- 
wise linear map with two control parameters: the slope and the bias. Their solution 
has two important features. First, it is complete, i.e. applicable for any (physically 
meaningful) values of the control parameters. Second, it reduces the problem of 
determining the transport coefficients to finding the sum of a quickly converging series. 
For the first time we thus have a model exhibiting irregular behaviour of the transport 
coefficients, for which a solution permitting an efficient numerical implementation is 
available. 

Klages and Klaufi [11] have recently used this exact solution to examine in detail 
irregular dependency of the drift velocity and diffusion coefficient on the control 
parameters of the system. Their aim was to verify an earlier hypothesis [12, 13] 
that the graphs of these quantities as functions of the map slope are so irregular that 
actually they form fractals. Using two numerical techniques: the box counting and the 
autocorrelation function methods, they found that the local fractal dimensions of these 
graphs are well-defined on small finite subintcrvals, but highly irregular functions of the 
control parameters. In other words they found that these graphs cannot be described 
with a single fractal dimension, but rather by a set of quickly varying local fractal 
dimensions. Taking this into account they put forward a hypothesis that the local 
fractal dimensions of these graphs as functions of the slope are fractal themselves. 

However, our recent numerical calculations [22] suggest that the Minkowski- 
Bouligand fractal dimension A of the graph of the diffusion coefficient as a function 
of the slope is equal 1, and that the convergence to this limit is slowed down by a 
logarithmic correction. Thanks to this logarithmic term the curve is a fractal. We 
also proposed a conjecture that the exponent controlling this correction depends on 
the slope of the map and equals either 1 or 2, depending on existence and detailed 
properties of a Markov partition. That A = 1 seems to be at odds with the above-cited 
results of [11], because the Minkowski-Bouligand dimension, if exists, is equivalent 
to the box-counting dimension [24]. However, one should note that while in our 
studies we focused on the point- wise dimension calculated only for those system control 
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parameters that correspond to finite Markov partitions, Klages and Klaufi computed 
the box-counting dimension on intervals of finite length, and these two values need 
not be the same. The aim of our present study is to support our earlier findings 
analytically. We also present new numerical results that confirm these results. 

The structure of the paper is as follows. Section 2 introduces the mathematical 
formalism. This includes the definition of the model, a brief description of 
the Minkowski-Bouligand dimension, the "oscillation" method of calculating it for 
continuous curves, and precise formulation of our hypothesis about the logarithmic 
corrections. Our main results arc contained in section 3, which presents a detailed 
calculation of the fractal dimension of the graph of the diffusion coefficient as a function 
of the slope for the case of zero bias. It also includes the derivation of the upper 
bound for the exponent controlling the logarithmic corrections. Section 4 contains a 
brief discussion of how our results can be generalized for the case of arbitrary bias and 
for other transport coefficients. Section 5 contains numerical results. These include 
analysis of the factors determining the value of the logarithmic correction as well as 
a study of self-similarity of the graphs. Finally, section 6 is devoted to discussion of 
results and conclusions. 

2. Basic formalism 

2.1. The model 

Following Groeneveld and Klages [10] we investigate a dynamical system with the 
equation of motion defined by a one-dimensional map M a ,b- R — > R parameterized by 
some real numbers a and b 

X n+1 = M a ,b(x n ) (1) 

where x n represents position of a particle and n is a discrete-time variable. The map 
is a piece-wise linear function given by 

M ^ = { a M a + b ( X + n)-n . ^R U) & 

where n is an arbitrary integer, and the slope a > 1 and the bias b £ R are the control 
parameters of the system (detailed analysis of the system symmetries [10] leads to the 
conclusion that the values of b can be restricted to — \ < b < |). The interval [—5,5) 
will be called "the fundamental interval" and denoted by Iq . However, instead of Iq 
one can choose Iq = (—5,5] as the fundamental interval in (2). It turns out that 
the transport coefficients (defined below) are the same whether we choose Iq or Iq 
as the fundamental interval. This equivalence has been employed in the derivation of 
the explicit forms of the transport coefficients [10]. 

The main features of the dynamical system defined by map M a ^ are depicted in 
figures la and lb. They show M a ^{x) for a = 2.83 and b = 0.1 and its action on two 
initial points xq. In particular, figure la presents the trajectory of xq = —\ using 
the fundamental interval [—5,5), while figure lb shows the trajectory of x = | with 
the fundamental interval (—5,^]. We chose these two trajectories because, as will 
be explained, they determine all transport properties of the system. The trajectories 
seem to be rather weakly correlated (which is typical for b ^ 0) and "random" . This 
randomness is related to the fact that, for nearly all values of a and 6, the iteration of 
the deterministic map M a ^ is equivalent to a stochastic Markov process of a random- 
walk type [17, 25, 26]: the deterministic trajectory of any point x n can be regarded 
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Figure 1. (a) The map M a ^(x) for particular values of the slope a = 2.83 and 
the bias b = 0.1 (thick solid lines); the dashed lines represent the first 33 terms 
of the trajectory of the point xq under subsequent iterations of M with xq = — i 
and the fundamental interval Iq . The arrow indicates the initial position xq . 
(b) The same as in (a), but for xq = ^ and the fundamental interval Iq . In both 
cases the transport coefficients are J fs 0.128 and D R3 0.0543. 



as a particular realisation of a random-walk process, and taking the average over the 
initial states Xq is equivalent to taking the averages over the corresponding Gibbs 
ensemble. For this reason the process defined by M a fc (or similar maps) is often called 
"deterministic diffusion". The two basic transport coefficients, the drift velocity J 
and diffusion constant D, are defined as 

j= lim M D= r lm < x l)-M\ (3) 

where (• • •} denotes the average over the uniform ensemble of initial values x . Figure 1 
illustrates also the role of using different fundamental intervals. Had we generated both 
trajectories using the same fundamental interval, they would differ only by a simple 
translation defined by (2). 

2. 2. Explicit form of the transport coefficients 

We will now briefly summarize Groeneveld and Klages's [10] method of calculating the 
transport coefficient for the model defined by the map M a ^. For each a > 1,6 G M, 
and e = ± one can define two infinite sequences of numbers, and n € r , r = 0, 1, . . ., 
consisting of reals and integers, respectively. Their values are uniquely determined by 
demanding that their first terms are given by 

I/o = |» «o = 0, (4) 
and for each r > 

K + Vr = ay%_ x +b, (5) 
with additional conditions 

n% e Z, and y% e Iq. (6) 
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Next we define "N-moments" : 



^k,i=N+-N^, (9) 

where fc, Z > 0. The basic transport coefficients, J and D can be now expressed ([10], 
cf. [23]) as 

j _ N 2fi _ N 3 ,q — N 2 .iJ + Ni. 2 J 2 . . 

In the particular case of b = we have, by symmetry, J = 0; hence the diffusion 
coefficient takes on a simpler form 

D=%&. (11) 

We will also use two important properties derived in Ref. [10]. First 

Ni,i > 0, (12) 

so that D and J are well-defined for any a ^ and b. Second, all numbers n% are 
bounded: 

K\<\a+\b\. (13) 

2. 3. Minkowski-Bouligand fractal dimension of the graph of a continuous function 

Let E be a bounded set in the plane. The e-Minkowski sausage of E, denoted by E(e), 
is the set of all the points whose distance to E is less then e. Let A(E(e)) be the area 
of E(e). Then the Minkowski-Bouligand fractal dimension is defined as [24] 

provided that this limit exists. In this case it is equivalent to the box-counting 
dimension [24]. A general relation between the Minkowski-Bouligand dimension A 
and the Hausdorff dimension Ajj reads A > Ajj. It is conjectured that for all strictly 
self-similar fractals A = A#. 

It has been recently proved [10] that the transport coefficients corresponding to 
the map M a ^ are continuous (but not diffcrentiable) functions of the slope a even 
though Nj t k are generally not continuous in a. Therefore we can calculate A for the 
graphs of D(a) or J(a) using the fact that the Minkowski-Bouligand dimension of a 
continuous real function f(f) defined on an interval [to,ti] can be evaluated through 
analysis of its Holder exponents [24] and r-oscillations [24] . The latter are defined as 

osc T (/;t)= sup f{t')- inf f(t'). (15) 

\t-t'\<T |t-*'l<r 

The Holder exponents are related to the oscillations as follows. If there exist constants 
c > and < H < 1 such that for all r 

osc T (/; t) < ct h (16) 
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then / is called a Holderian of exponent H at t. If constants c and H are independent 
of t then the fractal dimension A of / satisfies 

A<2-H. (17) 

Similarly, if there exist constants c > and H such that for all r 

osc T (f ;t) > ct h (18) 

then / is called an anti-Holdcrian of exponent H at t; if c and H are independent of 
t then 

A>2-H. (19) 

Relations (17) and (19) constitute a convenient means of calculating the Minkowski- 
Bouligand dimension A. Note that although they are supposed to hold "for all r", in 
practice it suffices to prove their validity only in the limit of r — > 0. 



2.4- Conjecture about the logarithmic correction to the fractal dimension 

In [22] we proposed the following conjecture: for the map M a _b in the limit of r — > 

osc T (D;a) [- log(r)] 7(a) , with C(o), 7 (a) > 0. (20) 

r 

Together with (15) - (19) this implies that the fractal dimension of the graph of 
the diffusion coefficient D as a function of the control parameter a is equal 1, but 
the convergence to this limiting value is logarithmically slow, with the exponent 7 
controlling the convergence rate. 

To justify this statement note that since lim r ^o(— log r) 7 /r ff_1 = 0forany7 > 
and H < 1, conjecture (20) implies that D(a) is Holderian of any H < 1. Using (17) 
we find that the fractal dimension of its graph is bounded from above by any number 
of the form 2 — H, with H < 1. Now we can take the limit of H | 1 to see that 
A cannot exceed 1. Similarly (20) implies that D(a) is anti-Holderian of exponent 
H = 1, and so A > 1. Hence A = 1. 

In the above analysis we have tacitly assumed that the coefficient C(a) in (20) 
can be bounded from below and above by positive constants independent of a. We 
will address this problem in section 5.2. 

Since D is a continuous, nowhere differentiable function of a [10], the fractal 
dimension A must be > 1 and exponent 7 must be strictly positive. This reasoning 
constitutes an alternative derivation of the lower bound for A and explains why 
7(a) > 0. Our numerical calculations (performed for the bias 6 = 0) indicated 
that exponent 7 depends on the slope a and is equal either to 1 or 2, depending 
on periodicity properties of the sequences y^ [22]. In particular, we conjectured that 
7 = 2 for the slopes a generating disjoint sequences y~£ and , and 7=1 otherwise. 

Note that even though A = 1, the graph of D(a) is not an ordinary smooth curve 
of dimension 1. Actually owing to the logarithmic term in (20), this graph is nowhere 
rectifiable (i.e., any of its arcs is "of infinite length"), which is typical of fractals [24]. 
Therefore the graph of D(a) locally resembles a well-known family of Takagi functions 
[6, 14, 15, 24], for their graphs are also nonrectifiable, their dimension is A = 1, and 
they satisfy (20) with 7(a) = 1 for all a. Takagi functions appear naturally in many 
contexts related to deterministic diffusion [6, 7, 14, 15]. 

Since hypothesis (20) opens a convenient way to explore fractal properties of 
graphs of the diffusion coefficient, its detailed analysis, carried out with both analytical 
and numerical methods, is one of the main goals of our present paper. 
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3. Fractal dimension of the diffusion coefficient for 6 = 

In the case of zero bias (6 = 0) the diffusion coefficient D vanishes for the slopes a < 2. 
Therefore we shall consider only the nontrivial case of 

a>2. (21) 

The basic properties of the diffusion coefficient are encoded in the sequences (n k ), 
introduced in Sec. 2.1. Before we will be able to tackle the problem of evaluating 
the fractal dimension of D(a), we need to derive several useful properties of these 
sequences. 

3.1. Basic properties of nk and yu for 6 = 

3.1.1. Symmetry. For 6 = the sequences n k and y k satisfy 

n k = vt = ~Vh . ( 22 ) 

so it will suffice to concentrate on the sequences n k and y k ; for simplicity we will 
denote their terms as n kl y kl respectively. These sequences may be either periodic or 
not. Periodic sequences nt (and yk) correspond to the so called Markov partitions of 
the interval 1^ (or 1^) [13, 14, 17, 25]. Parameters a for which sequences and yt 
are periodic will be called "Markov slopes". They form a dense set on (1, oo) [10]. 

3.1.2. Discontinuity in a. The numbers n^, as well as yk, can be considered as 
functions of a. Equation (5) implies that all functions rifc(a) and yk(a), fc = 1, 2, . . . , 
are discontinuous in a, which has a profound impact on the functional dependence 
of the transport coefficients J and D on a. This is illustrated in figure 2. Figure 2a 
presents the quantity of our primary interest: the diffusion coefficient D as a function 
of the slope a for 4 < a < 6 and for the vanishing bias 6 = 0. Although highly 
irregular, this function is continuous [10]. One might expect that so are -^3,0(0.) and 
TVi^a), because their quotient equals to D. However, as depicted in figures 2b and 2c, 
actually they are discontinuous. From (7) - (9) we conclude that this must be related 
to discontinuity of functions n k (a). Two of them, r^(a) and n^(a), are shown in 
figures 2d and 2c, respectively (as horizontal line segments). For completeness, in the 
same figures are also depicted the graphs of y^~(a) and y^(a). As expected, there 
is a clear correlation between discontinuities of rifc(a), yk(a) and discontinuities of 
-^3,o(a), A^i ; i(a). Moreover, there is also a correlation between these discontinuities 
and "irregularities" of D(a). 

The union of all discontinuity points of nk{a) for all k will be denoted by 5 and 
we will introduce H as the complement of S on the interval (l,oo). Importance of 
these two sets is related to the fact that functions Nk,i(a) are continuous on S and 
discontinuous on S. Moreover, functions (a) and yk(a) are also continuous on S 
(however, in contrast to Nk,i(a), they can be also continuous on some points of S). 
Both S and S are dense on (1, oo) [10]. 

3.1.3. Approximate linearity of yk(a). Figures 2d and 2e suggest that 2/2(0) and 
2/3 (a) consist of almost linear segments even though in reality they are polynomials 
of degree 2 and 3, respectively. Moreover, the density of discontinuity points of n 3 
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Figure 2. Graphs of the basic quantities used in the mathematical derivation as 
functions of the slope a, calculated for the bias 6 = 0. (a) The diffusion coefficient 
D = N:ifi/Ni \; (b) -/V30; (c) Nn; (d) ri2 (horizontal, thick lines) and ?/2 (thinner 
curves); (e) nz (horizontal, thick lines) and 1/3 (thinner, almost linear curves). 
Except for D, all of them are discontinuous. Nj^ are continuous on a dense set S 
and discontinuous on a dense set S which is the union of all discontinuity points 
of ni(a), n2(a), no,{a) . . . (see text for details). 



is considerably greater than that of n2- These two important observations can be 
justified in a more formal way. First of all note that (5) implies 

fe-i 

y'k = Vk-i + ay'k-i = Vk-1-ja 3 , (23) 

3=0 
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where we used a short-hand notation y' k = dy k (a)/da. Therefore, since |j/fe(a)| < \ 
and yo = \, the derivative dy k (a)/da is bounded 

a k - 1 (a-2) + l . a k -l 

2(a-l) * * * 2(a-Ty (24) 
Taking into account (21) we conclude that for large k 

J/fc-a*" 1 . (25) 
Since the range of each y k is limited to [— |, |], the distance between two consecutive 
discontinuity points of yk(a), denoted as 6 k (a), is of order < a 1_fe , i.e. decreases rapidly 
with k. Similar reasoning leads to the conclusion that the second derivative of yk(a) is 
of order < 2(fc— l)a . Thus the maximal error resulting from truncating the Taylor 
series of yk(a) after the linear term for the displacement 6k (a) ~ a 1 ~ fe is of order 
(k — l)/a k ~ 1 and vanishes quickly as k goes to infinity. Therefore functions yk{o) can 
actually be approximated to a high accuracy by piece-wise linear segments. 

3.1.4- The limits of a going to ao from below and from above. Although functions 
Nk,i(a) are discontinuous on S, they are left continuous. Moreover, it turns out that 
for any a £ (l,oo) the limit of Nkj(a) for a going to a from above is also well- 
defined. We will denote it by Nk,i(a>o)- Its value can be related to Nk,i(ao) as follows. 
For ao <E H, by definition of S, Nk,l(ao) = N k j(ao). On the other hand, for a £ 5 
the sequence no(ao), ni(ao), ... is strictly periodic. Denoting the period length by L 
we thus have nk (ao) = nk+L{ao) for all k > 0. For a tending to ao from above the 
sequence rife (a) will converge to a new sequence hk(ao) 

' nfe(ao), k < L 

nk(a ) = < n k (a a ) + 1, k = L (26) 
_ -nfc(ao), fc > L 

Now Nk,l(do) = li m o-»o + Nk,i{a ) can be determined upon substituting n^(ao) for n k 
in (7) - (9). Functions Nk,i(a) are right continuous. As argued in [10], all transport 
coefficients of the system defined by the map M a ,b are continuous in a. Consequently, 
functions Nk,i(ao) can replace Nk,i(a) in (10) to calculate the transport coefficients of 
the system. It is also interesting to note that even though the functions Nk,i(a) and 
Nk,i{oi) are different, their graphs must be the same. For this reason figure 2b actually 
shows both Ns t o(a) and A r 3 j o(a). Similar conclusion pertains also to figure 2c. 

3.2. Oscillations at elements ofE 

Let ao be an arbitrary element of S. This means that all functions yk{o) are 
continuous at ao, and so for each integer k there exists e > such that yo, j/i, • • • , yk 
and no,ni, . . . ,nk are continuous and diffcrentiable functions of a on the interval 
(ao — £o, ao+£o). Therefore for any a G (ao — £o, ao+£o) all functions nj(a) up to j = k 
are independent of a. Let < r < e . Since yk(a) is continuous on (a — £o, a + £o), 
we can write down its Taylor expansion about ao and, following our discussion in 
section 3.1.3, truncate it after the linear term, obtaining yfe(ao + r) w yk{o-o) +i~y' k (ao). 

Let fj,(r; a) be the maximum integer j such that yj(a) and Uj(a) are continuous 
on (a — r,a + r) . This means that for any ao £ S and a fixed \i numbers \i and 
r are related through y M (a ) + ry^(a ) « \ or y M (a ) - ry^(a ) w -| (whichever 
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gives smaller r). This, together with (25), leads to the following implicit equation for 

M (r;a)«log a i-fe^. (27) 

There are two cases: either the sequence yo(do), J/i(ao), ■• ■ is periodic (i.e. the 
system has a Markov partition) or not. In the former case the numerator in (27) is 
bounded from above and below by positive numbers; in the limit of r — > we thus 
have 

/u(r; a) « log a t _1 + const (for a <G 5). (28) 

Clearly /i(t; a) — > oo as t — > 0. In other words, in the limit t — > (which is essential 
for the analysis of Holder exponents), the number of functions in the sequences yk{a) 
and nk (a) that are continuous on (a — r, a + r) goes to infinity. 

The case where the sequence yo(ao), yi(ao), . . . is not periodic is far more 
complicated. The main problem comes from the fact that now the numerator in (27) 
can get arbitrarily close to zero, which may invalidate (28). In the present approach 
we restrict ourselves to the cases where the system has a Markov partition and so the 
sequence yo{ao),y\(a,o), ■ ■ ■ is periodic. 

As expressed by (11), the diffusion coefficient is the quotient of two quantities, 
A3.0 and Ni t i. However, by a straightforward calculation we can verify that for any 
t > and any functions f(t) and g(t) bounded on [t — T,t + r] (no matter if continuous 
or not, but the lower bound for \g\ must be ^ 0) 

osc T (//s) = max T (//5) - mm T (f/g) 

osc r (/) max r (g) + osc T (g) min T (/) 
min r (g) max T (#) 
< CiOSC r (/) + c 2 osc T (5) (29) 
where we used a short-hand notation max T (/) = supu_ t /i <r f(t') and min r (/) = 
inf \t-t'\<r f(t'), while c\ and c 2 are some nonnegative parameters independent of r. 

Thanks to (12) we can apply (29) to D(a), substituting A^o for /, Ni.i for g, 
and a for t. The problem of finding the r-oscillations of the diffusion coefficient D can 
be thus reduced to that of finding osc^A^o; a) and osc r (Ai j i; a). 

Equation (8) implies that for arbitrary integers k, I and a = ao ± r 

N kil (a )-N kil (a)<^ {a r [N+(a )] k a-[A+(a)] fe } r\ (30) 

r=/i(T;a ) 

with \x related to r through equation (28). Here we have employed the fact that, by 
definition of /U, n r (a) = n r (a ) and hence A+(a) = A+(ao) for all r < fi(r;a). Using 
(13) we conclude that 

\N+\<ctr + c* 2 , (31) 

with c\ > 0, c* 2 > being some constants. It will be convenient to assume that c\ 
denotes the smallest possible value satisfying (31), c\ = limsup,,^^ \N+(ao)\/r. 

Thanks to (31) an upper (as well as lower) bound for Nk,i(ao) — Nkj(a) can be 
found by substituting arithmetic sequences for A+(a ) and A+(a) in (30). Their 
common differences may be different, but the terms for r = [i must be the same or 
differ by a small number less than + cf (13). This leads to 

\Nk,i(a) - N k ,i(ao)\ < a -^ k+l -\ for c\ > 0, (32) 

\N k ,i(a)- N kil (ao)\ < a"V, for c\ = 0. (33) 
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This important result can be also obtained by noticing that the major contribution 
to the sum in (30) comes from its first few terms. Hence, using (28), 

|JV fc>I (ao±r)-JV fc ,,(ao)| < rln^" 1 ^" 1 ), c\ > 0, (34) 
\N kil (a ±T)-N k j(a )\<Tlii l (T- 1 ), c\ = 0, (35) 

which implies 

osc r (7V M ;a) < t ln*^ 1 (r^ 1 ), c\ > 0, (36) 
osc T (N kX ,a)<Tln l (T- 1 ), c{ = 0. (37) 

Consequently, by (29) and (11), 

osc r (L>;a) < Tin 2 ^ 1 ), c{ > 0, (38) 
osc r (^;a) ^lnfr 1 ), c?=0. (39) 

for all a £ H. 

Equations (38) and (39) give the upper bound for the exponent 7 controlling the 
logarithmic correction to the local fractal dimension. As we can see, this upper bound 
depends on c\. This constant has a simple physical interpretation. If the trajectories 
of points xo = ±i under M a ^ (like those depicted in figure 1) remain confined inside 
a neighborhood of the origin, we can take c\ = 0, which implies < 7 < 1. Otherwise 
c\ > and < 7 < 2. 



3.3. Oscillations at elements o/S 

The argumentation presented in section 3.2 might seem of little use for elements of 
the set S, at which functions N ky i(a) are discontinuous. However functions N k d(a) 
are left continuous and have a well-defined limit from the right, Nkj(a). Functions 
Nk,i(a) are right continuous and, as argued in section 3.1.3, can replace Nk,i(a) in 
calculations of the transport coefficients. This suggests a method of circumventing 
the problems related to discontinuity. The oscillation of the diffusion coefficient D on 
the interval [ao — r, ao+r] with ao <E S can be calculated as the maximum of two values: 
the oscillation of (A^o/A^i) (a) on [a — r, a] and the oscillation of (iV^o/iV^i) (a) on 
[a, a + r]. In each case we can now apply the reasoning of section 3.2 (its most delicate 
point, validity of approximation (28), can be justified rigorously because the sequence 
rife (a) is periodic for all a £ S). This implies that (38) and (39) are valid also for all 
a £ S, which suggests that the fractal dimension of the diffusion coefficient is 1 for all 
Markov slopes a £ (1, 00). 



4. The case b^0 

Besides the fractal properties of the diffusion coefficient D, in the asymmetric case 
i) / we can investigate also those of the drift J, which for b ^ is also a highly 
irregular function of the control parameters [10]. It turns out that the methods 
developed in section 3 can be readily extended to the case of nonvanishing bias b 
and to the analysis of arbitrary transport coefficients. Therefore we will only discuss 
the major differences and difficulties appearing when applying these methods to the 
more general case 6^0. 

First, although derivation of the basic relation (27) remains practically 
unchanged, we see that restricting our calculations to the case 1 < a < 2 is fully 
acceptable only for b = 0. This is related to the fact that for b ^ the transport 
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coefficients can assume nontrivial values also for a < 2 [10], so this parameter region 
also deserves attention. The assumption a > 2 greatly simplifies the proof of (25) 
(cf. (24)); without it we have to take into account correlations between consecutive 
terms of the sequences . Our numerical analysis shows that (24), and hence (27), 
is valid also for 1 < a < 2 (at least in the limit of k — > oo), but we failed to find a 
rigorous proof of it. 

Second, the sequences and are no longer related to each other by (22) and 
in practice must be considered as independent of each other. This, fortunately, is not 
a serious difficulty - actually relation (22) was used mainly to simplify notation. 

Third, just as in the case b = 0, we can justify transition from (27) to (28) only 
in the case where both n\ and are periodic, i.e. when the system has a Markov 
partition. 

Fourth, we still need an analytical argument justifying our claim that the graph 
of (D;a) is uniformly anti-Holderian for all b. 

Although in section 3 we focused on the fractal properties of the diffusion 
coefficient D, our analysis was as general as possible. In particular inequalities (36) 
and (37) were derived for functions Nk.i(a) with arbitrary k,l. Since all transport 
coefficients, including the drift velocity J, can be expressed as functions of the 
quotients Nk,i /N\ t \ [10] , the methods of section 3 can be applied for arbitrary transport 
coefficients. 

Consequently, we can follow the ideas of section 3 to argue that the local fractal 
dimension of the graphs of transport coefficients as functions of a is 1 for the control 
parameters a and b e R such that a > 2 and the system has a Markov partition. 
Based on our numerical calculations [22] and heuristic arguments we conjecture the 
same fractal behaviour for all other pairs (a, b). 

5. Numerical results 

5.1. The role of c\ 

Equations (38) and (39) imply that the exponent 7 controlling the logarithmic 
correction is bounded from above by 2 for c\ > and by 1 for c\ = 0. On the other 
hand, however, in [22] we conjectured that the value of 7 depends on properties of 
the periods of sequences y~^ and y^ (i.e., whether they consist of the same or different 
terms). A question arises whether the two approaches lead to the same conclusions 
for 7 and if not, what does its value really depend on? 

To answer it we will study a particular case of the bias 6 = and the slope a 
equal to the largest root of a 5 — 2a 4 — 2a 3 + 2a 2 + 1, i.e. a « 2.455. The trajectory 
of Xq = \ is xf w 1.23, x\ rts 1.56, rts 0.92, x\ w 0.80, x\ = Xq, . . ., and so the 
sequence x+ is periodic, x^ +m = x+ . The trajectory of Xq is given by x~ n = — x+ , 
and hence is periodic, too. Since all terms of the sequence xf, are different from those 
of , according to [22] the value of 7 should be 2. This is in conflict with (39), as 
the trajectories of x = ±5 are bounded, hence c\ — and, following (39), 7 < 1. We 
checked this numerically using the arbitrary precision arithmetic library GMP [27]. 
Our results are presented in figure 3. It clearly suggests that 7 = 1. This implies 
that our earlier hypothesis about the logarithmic correction exponent 7 was incorrect. 
Our present, deeper insight into mathematical nuances of the problem leads to the 
conclusion that the main factor determining the value of 7 (at least for systems with 
periodic trajectories of xq = ±|) is whether the constant c\ vanishes or not. 
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Figure 3. Oscillations of the diffusion coefficient D for b = and a equal to the 
largest root of a 5 — 2a 4 — 2a 3 + 2a 2 + 1, i.e. a »J 2.455, rescaled in accordance 
with (20) and calculated on intervals (a-r,a + r) for 10" 100 < r < 10~ 2 . The 
inset presents a blowup of the results and reflects periodicity of the curve. 



Figure 3, and especially the inset in it, indicates also that in the case of the map 
M aj {, oscillations of the diffusion coefficient really deserve their name. Apparently the 
graph of osc T (D)/r as a function of logr -1 is a composition of a linear and periodic, 
"oscillating" function. The period length A of these "oscillations of oscillations" can 
be estimated using (28), 

X = a L , (40) 

where L is the period of y^\ in our case L = 5, A = a 5 « 89, log A « 1.95, in excellent 
agreement with the data in figure 3. 

5.2. The uniformly Holder and anti- Holder conditions for D(a) 

As we explained in section 2.4, analysis of Holder exponents can be useful in 
exploration of fractal properties of a curve if this curve is uniformly Holderian (or 
anti-Holdcrian) , i.e. if the coefficient C(a) in (20) can be bounded from below (or 
above) by a constant number. We decided to check these properties numerically for the 
vanishing bias b, as in this particular case all Markov slopes a are algebraic numbers. 
They can be thus generated in a systematic way by first generating polynomials 
with integer coefficients and then identifying Markov slopes with the largest roots 
of the polynomials. Moreover, the coefficients of each such polynomial determine the 
numbers n%, hence whether the constant c\ vanishes or not, and thus whether the 
exponent 7(a) in our conjecture (20) is expected to equal 1 or 2, respectively. 

To check whether D(a) is uniformly anti-Holderian we concentrated on the case 
7(a) = 1, i.e. c\ = 0. Using a computer program we found all Markov slopes a € (2; 8] 
that are algebraic numbers of degree d < 4 and correspond to c\ = 0. Then, for each 
such a slope, we calculated the coefficient C{a). To this end we assumed that at each 
investigated slope the graph of osc T (D; a) as a function of log t" 1 behaves like the one 
depicted in figure 3. To increase accuracy we eliminated local periodic oscillations of 
the curves by using (40). As expected, the points of each curve at abscissas forming 
an arithmetic sequence with a common difference log A turned out to form practically 
a straight line with the slope quickly converging; we identified this limiting slope 
with C(a). In each case the sum of squares of the residuals from the best-fit line 
was less than 10~ 10 . The results obtained in this way are depicted in figure 4a. As 
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Figure 4. (a) The coefficient C(a) for all Markov slopes a £ (2; 8] that correspond 
to 7 = 1 and are algebraic numbers of degree d < 4 (2054 data points). The inset 
presents the same quantity calculated for a S (2; 3] and d < 6 (439 data points). 
The plus (+) and cross (x) signs mark the estimated values of C(a) for non- 
Markov slopes a = it and a = 4.5, respectively, (b) The same as in (a), but 
for 7 = 2 and d = 1 (stars), d = 2 (crosses), and (i = 3 (small squares). The 
logarithmic scale for the ordinate axis was used to visualize the structure of C(a). 



expected, the density of the points increases with a. Moreover, apparently they lie on 
a continuous curve that crosses the x-axis only at a = 2. These findings arc supported 
by the results depicted in the inset of this figure and obtained for all Markov slopes 
a E (2; 3] that correspond to c\ — and d < 6. That C(a) tends to a continuous curve 
implies that the graph of D as a function of a is uniformly anti-Holderian with H = 1 
and 7 = 1 at least on the set of all Markov slopes a > 2 for which cf = 0. Note that 
this set is dense on (2, oo). 

To check whether D(a) is uniformly Holderian we focused on the Markov slopes 
with 7(a) = 2, i.e., c\ ^ 0. We investigated numerically all Markov slopes a G (2; 8] 
that are algebraic numbers of degree d < 3 and correspond to c\ ^ (286 data points). 
We found that at each such a slope the exponent 7 is actually equal to 2. However, as 
seen in figure 4b, the coefficient C(a) turned out to be a highly irregular, discontinuous 
function of a. Fortunately, this function is apparently bounded from above. Taking 
into account our earlier result for c\ = 0, we conclude that the graph of the diffusion 
coefficient D as a function of the slope a is uniformly Holderian (with H — 1 and 
7 = 2) at least on the (dense) set of all Markov slopes. 

The analysis of the Holder and anti- Holder condition at non-Markov slopes is far 
more difficult. Such values of a were already investigated numerically in [22]. The 
conclusion was that in this case most probably 7=1, but this property cannot be 
definitely established numerically due to very poor convergence. In several cases there 
is a clear linear trend in graphs of osc r (D; a) as a function of logT -1 , which permits 
to estimate the value of C(a) (assuming 7(a) = 1). For example, using the data of [22], 
we can estimate C(tt) « 0.215 and C(4.5) w 0.15. As seen in figure 4, these values are 
very close to those obtained for the nearby Markov slopes. However, for many other 
values of a no reasonable numerical estimation of C(a) is possible. 



5.3. Oscillations of Nk,i(a) 



Results of section 3.2 reveal that the logarithmic corrections to fractal dimensions of 
graphs of transport coefficients must be attributed to similar logarithmic corrections 
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Figure 5. Oscillations of JVi i and ./V3.0 rescaled in accordance with (41) 
calculated on intervals (a — r, a+r) for 10 — 100 < r < 10~ 2 , with 6 = and a equal 
to the largest root of (a) a 5 - 4a 4 + 2a 3 + 3a 2 - 2a + 4, (b) a 5 - 4a 4 - 4a 2 + a + 4. 
Results for A^q in the right panel were fitted by a quadratic (a solid line, hardly 
distinguishable from numerical data). The insets are blowups of the curves to 
show their periodicity and discontinuity in a. 

in functions N k ^(a). For a £ 5 these latter functions are continuous in a and we 
expect them to satisfy (20) with Nk,i substituted for D, i.e., for r — > 

sc r (Ar fc , ; ;a)^ [log(r _ 1)]7M(o ) (4l) 

where, following (36) and (37), 

< 7fc,j < k + 1 - 1 for cj > 0, (42) 

< 7 fei ; < / for = 0. (43) 

This conjecture relates the fractal nature of transport coefficients to fractality of 
functions Nk,i(a) which — although of no immediate physical meaning — have much 
simpler form and are more amenable to rigorous treatment. 

To verify (41) - (43) we have calculated oscillations of A^o and N\ t i for the 
vanishing bias b and two values of the slope a: the largest root of a 5 — 4a 4 + 2a 3 + 
3a 2 - 2a + 4 (i.e., for a w 3.04) and the largest root of a 5 - 4a 4 - 4a 2 + a + 4 (i.e., 
for a w 4.20). For these parameter values the system has a Markov partition [22]; 
moreover, the constant c\ — for a w 3.04 and c\ > for a w 4.20. 

Our results for a w 3.04 and a w 4.20 are depicted in figures 5a and 5b, 
respectively. As we can see, our data clearly indicate that in both cases 71,1 = 1. 
However, 7 3j0 = in figure 5a (i.e., for c* = 0) and 73^ = 2 in figure 5b (i.e., for 
c\ > 0). All these results are in perfect agreement with (43) and actually correspond 
to the upper bounds implied by these equations. Note that for c\ = equation (43) 
actually gives not only the upper bound, but the exact value of 7/^0 = for all k > 0. 

The insets in figure 5 suggest that, just as in the case of graphs of osc T (Z?)/r 
discussed in section 5.1, on finer scales the graphs of osc T (N k j)/r as functions of 
logr -1 have periodic components. Their period lengths A can be calculated using 
(28), which gives A = a 6 w 794 for a w 3.04 (figure 5a) and A = a 4 w 311 for a w 4.20 
(figure 5b). These theoretical values are in excellent agreement with numerical data 
used to produce figure 5. 
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Figure 6. (a) The graph of the diffusion coefficient obtained for the slope ao fs 
2.455 and the bias 6 = (dotted line) and its consecutive A-fold magnifications 
about the point (a ,D ), where A = « 89 and D = D(a ) » 0.0967. (b) The 
same as in (a), except that now each A-fold magnification has been followed by 
a linear transformation (x,y) i— > (x,y — sx) with s ra —0.1969; all 4 curves have 
been plotted, but those obtained for m = 1,2,3 differ from each other by less 
than the line width. 



5.4. Self- similarity of the graphs 

In the limit of r — > the periodicity of the curves presented in figures 3 and 5 becomes 
practically perfect. This phenomenon could be a consequence of local self-similarity 
of functions D(a) and Nk.i(a) under enlargement of their graphs by a scale factor 
A. To verify this hypothesis in figure 6a we have plotted the graph of the diffusion 
coefficient D as a function of the slope a, with the origin moved to the point ao ~ 2.455 
and Dq = D(ao) ~ 0.0967. Then we plotted three consecutive blowups of this graph 
obtained by magnifying it about the new origin by scale factors A, A 2 , A 3 , where 
A = Oq « 89. As we can see in figure 6a, this simple transformation of graphs does 
not lead to perfectly self-similar structures; nevertheless, the curves thus obtained do 
exhibit some kind of similarity. 

The effect of actual self-similarity can be achieved if each A-fold enlargement will 
be followed by a correction transformation related to the nonvanishing logarithmic 
correction exponent 7. The 2-step transformation T can be expressed as 

T(x,y) = (Xx^y-miXx)), (44) 

where A is the scale factor and K\ is a correction term. Our numerical results 
suggest that for 7 = 1 this complementary step is a simple linear transformation, 
K\(x) — sx, with s being a constant depending on the slope a, the bias b and the 
function the self-similarity transformation T is applied to. In the case of the data 
shown in figure 6 the constant s w —0.1969. The net effect of this two-step scaling 
procedure applied iteratively m — 0,1,2,3 times is depicted in figure 6b. As we 
can see, the convergence is very fast and the curves obtained after m = 1, 2 and 3 
elementary 2-step transformations are already practically indistinguishable from each 
other. Closer analysis of the numerical data for m = 0, . . . , 10 revealed that with 
each 2-step transformation we get closer to the self-similar limit curve by about two 
significant figures. Similar behaviour was found for several other points a , both for 
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Figure 7. (a) Correction functions rei and K2 used in self-similarity 
transformation (46) for the graph of the diffusion coefficient D(a) about the point 
(ao,Do) with ao ~ 4.20, Do = D(ao), and the bias b = 0. (b) The self-similarity 
effect of applying T m , as described by (46), to the graph of D(a), for the same 
parameter values as those used to generate panel (a). Note that the dotted line 
(corresponding to m = 0) actually represents the graph of the original function 
D(a) translated by a vector (— ao, —-Do) and x can be interpreted as A m (a — ao). 
The graphs obtained for m > 1 are practically indistinguishable from each other. 

graphs of the diffusion coefficient D(a) and the auxiliary functions N\ t i(a), N^fi(a), 
provided that the corresponding logarithmic correction exponent 7 was 1. Note that 
transformation of type (44) holds also for Takagi functions. 

For 7 = 1 the net effect of applying consecutively m transformation T is 

T m (x,y) = (\ m x,\ m y-m Kl (X m x)), (45) 

where linearity in m of the correction term corresponds to 7 = 1. For 7 = 2 the 
situation gets more complicated. We now expect that the correction term after m 
self-similarity transformations is quadratic in m, 

T m (x, y) = {\ m x, \ m y - m 2 K 2 {\ m x) - m Kl (X m x)), (46) 

with some real functions k\ and n 2 . A side-effect of such a form of the correction term 
is that, in contrast to the case 7=1, now the second step in each elementary 2-step 
transformation depends on m and hence slightly differs from each other. We have 
confirmed (46) numerically, finding the convergence to be very fast. In particular, 
figure 7 presents results obtained for the slope ao « 4.20 and bias 60 = (these 
parameters were already used in figure 5b). On the left panel we show correction 
functions K\ and k 2 . As we can see, k 2 is linear, while K\ is an irregular (perhaps 
fractal) function of its argument. Although the graph of K\{x) in figure 7a is quite 
similar to that of D(a) (dotted line in figure 7b), these two functions are different. 
On the right panel we present effects of applying T m , m = 0, 1, 2,3, to the graph of 
D(a) in the vicinity of the point (a , D(a )). The convergence to the self-similar limit 
is very fast. Actually for m > 1 the graphs of T m in figure 7b differ from each other 
by less than the line width. 
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6. Summary and conclusions 

We have analysed both analytically and numerically the local (Minkowski-Bouligand) 
fractal dimension of graphs of transport coefficients as functions of the slope a in a 
simple dynamical system defined by a piece- wise linear map M a ^. We showed that for 
all values of the control parameters a and b such that the system defined by M a ^ has 
a Markov partition, the dimension of such graphs is 1, probably with a logarithmic 
correction. This correction renders the curves nowhere rectifiable (i.e., any of their 
arcs is "of infinite length"), which is typical of fractals. Our results imply also that 
the Hausdorff dimension of these curves is 1, too. 

We have also found the upper bound for the exponent 7 controlling the logarithmic 
correction. This bound turned out to be either 1 or 2, depending on whether a 
constant c\ vanishes or not, or, equivalently, whether the trajectories of xq = ±| 
remain bounded or not. This finding extends our earlier conjecture [22] about factors 
determining the value of 7. 

Using numerical calculations we have shown that at Markov slopes a the graph 
of the diffusion coefficient is actually self-similar under the action of a special, two- 
step scaling transformation. The first step is a simple uniform enlargement by a scale 
factor A. This transformation must be followed by the second step, a "correction" 
transformation of the graph, which is related to the non-zero value of the exponent 
7 controlling the logarithmic correction to the local fractal dimension. We found the 
precise formula for the value of the scale factor A as well as a general form of the 
correction transformations. In particular, we found that the scale factor A depends on 
the slope a and the periodicity of the trajectories of points x = ±i. The correction 
transformations turned out linear (affinc) for 7 = 1 and different, more complicated 
for 7 = 2. At non-Markov slopes the graphs are probably not self-similar. 

The graphs of transport coefficient have thus surprisingly reach spectrum of local 
fractal properties. Our numerical results suggest that the curves are self-similar on a 
dense set of points (at which the system has a Markov partition), but at each such 
a point the scale factor and the correction transformation are different. In this sense 
the curves are multifractals. 

Our analytical approach has, however, some limitations. The most serious one 
is related to the fact that it is applicable only to systems in which the trajectories of 
points xq = ±5 are periodic (or, equivalently, to systems with a Markov partition). 
We believe that the typical value of 7 at non-Markovian slopes a is 1, but in general 
this value can take on any value between 1 and 2. This conjecture, however, requires 
further studies. Proving that the investigated curves are uniformly Holderian with 
H = 1 and 7 = 2 and anti-Holdcrian with H — 1 and 7 = 1 is another open problem. 
So far we have been able to justify these properties only numerically. 

Numerical results suggest that our analytical approach gives not only the upper 
bound, but the exact value of 7. To prove this hypothesis rigorously would require 
finding a "tight" lower bound for 7. However, as is well known [24], finding the lower 
bound of a fractal dimension is usually much more difficult that finding the upper one. 
This is also true in our case. For example, when we were looking for the upper bound 
of oscillations of the quotient of two functions, we did not have to take into account the 
fact that actually their values were correlated. Applying the same general assumption 
for the lower bound of 7 yields only a trivial result 7 > 0. Therefore finding a better 
lower bound for the logarithmic correction is another problem for future studies. 

Another question is how to reconcile our results with those obtained by Klagcs 
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and Klaufi [11], who analysed the box-counting dimension of D(a) on intervals of finite 
length and found it to be greater than 1. One possible way of answering this problem 
could be to analyse how the local self-similarity of the curves depends on the box size. 

Although most of our analysis was carried out explicitly for the graph of the 
diffusion coefficient D as a function of the slope a for a vanishing bias 6, we showed 
how it can be generalized for other transport coefficients and for arbitrary bias b. 
Similarly we expect that the main conclusions derived here for the simple piece- 
wise linear map M a .b remain valid also in the case of more realistic (and complex) 
dynamical systems. After all, in our model the main reason for irregular, fractal 
dependence of the transport coefficients on the slope a is an extreme sensitivity of the 
(periodic) Markov orbits on the control parameters; similar sensitivity is observed in 
more complex models, e.g. the multiBakcr or the Lorentz gas [1, 2, 3, 6, 17]. Transport 
coefficients in these models were already called "fractal", but only in the sense of 
"unexpectedly irregular". We hope that our present approach will help justifying a 
conjecture that the transport coefficients in those more realistic models form fractals 
in a strictly mathematical sense of this term. 
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